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Abstract 



Explicit formulae are given for the saddle connection for an inte- 
grable family of standard maps studied by Suris. A generalization of 
Melnikov's method shows that, upon perturbation, this connection is 
destroyed. We give explicit formula for the first order approximation 
of the area of the lobes of the resultant turnstile. It is shown that 
the lobe area is exponentially small in the limit when the Suris map 
approaches the trivial twist map. 

AMS classification scheme numbers: 

34C35,34C37,58F05,70H15,70K99 
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1 Introduction 

Standard maps are area-preserving diffeomorphisms of T x M given by 

f^e,r) = {e + r + v'{e),r + v'{e)), (i) 

where the potential, V, is periodic, V{6 + 1) = V{6). The case where 

F(e) = ^cos(27re) (2) 

is known as the standard or Taylor-Chirikov map. This model is important 
because it gives a local description of nonintegrable two degree of freedom 
Hamiltonian dynamics. 

Twist maps, of which the standard map is an example, will be our major 
concern in this paper (for review, see [17]). Such maps have Lagrangian 
generating functions, S{6,6'), which generate the map implicitly through 
the equations 

r= -diS{e,e'), 

r'= d2S{e,e'). (3) 



To generate a map, S must satisfy the twist condition that the second equa- 
tion above can be inverted to obtain 6'{r, 6); this occurs, e.g. if 

81828 < , (4) 

and implies the geometric condition that vertical lines tilt to the right upon 
iteration, ^ > 0. Furthermore, we assume that our twist map has zero net 
flux, which is equivalent to 

s{e + i,9' + 1) = sie,e') . (5) 

For maps of the standard form, the generating function is 

s{e,e') = ^{d'-df + v{e). (6) 

In this paper we discuss a standard map introduced by Suris [21, 17] 
with the potential given by 

Vie)= -|/o^rfitan-i(T5i^I&)) 

= ^SR[dilog(l + 5)-dilog(l + 5e2^^^)] , 
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Figure 1: Some contours of Z^, for 5 = 1/3. 

where the dilogarithm is defined by 

diiog(^) ^ r ^^dz . 

Ji I — z 

The Suris map, Z^, i.e., the standard map with this potential, is inte- 
grable with integral 

f) = COS 7rr + 5 cos 7r(20 — r) , 

i.e. Ig o fg = Ig. Contours of Ig are shown in figure 1. The map is integrable 
for any 5, however, we will consider the case < ^ < 1, as the topology of 
the saddle connections changes at 6 = 1. For < 6 < 1 the map fg has 
hyperbolic fixed points at Za = (—^,0) and z^ = (^,0), that are connected 
by two saddle connections, forming the upper and lower separatrices of the 
fixed point resonance. 

Since the Suris map is integrable and has twist, Birkhoff's theorem [17] 
implies that the saddle connection between the equivalent points Za and Zf,, 
which is a rotational invariant circle, is necessarily the graph of a function 
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r = x(^)) for 6 G = [— |, |]. The dynamics on the saddle connection from 
(1), gives the circle diffeomorphism 9 h{9) = 6 + x{G) + V'{9). 

We will show that, under perturbation, the stable and unstable manifolds 
of the perturbed twist map intersect transversally. We use a modification of 
Melnikov's method [18] for twist maps [7, 8, 9]. 

We perturb the system by adding to the original potential any peri- 
odic function P{9). We assume for simplicity that P also satisfies P(— |) = 
P(i) = and P'{-^) = P'il) = 0. Then for small enough e, Ss + eP is 
a twist generating function, and the corresponding map fg^e has hyperbolic 
fixed points at Za and zi,. As was shown in [8], the Melnikov series, given by 



converges absolutely and imiformly to a function on compact subsets of 
U. Furthermore, if the function L has a nondegenerate critical point in U, 
then the manifolds W'"{za, fs,t) and W^{zh, fs,e) intersect transversally for e 
small enough (the same conclusion is valid for W'^{za, fd,^) and W'^{zi,, fd,^))- 

Thus the function L provides sufficient conditions for the transversal 
destruction of the saddle connection, just as in the classic applications of 
the Melnikov integral. In section 2 we formulate a slightly stronger version 
of the previous result (theorem 2). 

Below we consider the case P{6) = cos^ tt6. We show that for e > and 
small enough there are two distinct transversal heteroclinic orbits from Za to 
zi,. These delineate a turnstile with two lobes that define the areas that are 
transported into and out of the fixed point resonance upon each iteration of 
the map. The area of the lobe gives a coordinate independent measure of 
the separation of the perturbed stable and unstable manifolds as well as the 
flux from one region to another. 

The exact, order e approximation for the lobe area is obtained in section 4 
and compared with numerical results in section 7. We summarize our results, 
as 

Theorem 1 (Main Theorem) Let 



oo 



2 

Vs{6) = / arctan 

vr Jo 



1 + Jcos(27rt) 



5 sm{2nt) 



dt 



(7) 
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be a perturbed generating function with corresponding twist map fg^^. Then 
for all < 6 < 1 and e small enough there are two orbits heterodinic from 
Za = (—^,0) to Zb = (^,0). The lobe defined by these two orbits has an area 
A given by 

A{S,e) =eT{u) + Oie^) (8) 

where 

rM-i + 8i:UI^, (9) 

k=l 

and u = {l- \/S)/{l + ^^). 

Note that is a natural parameter for the stable and unstable manifolds, 
and hence for the area of the lobe, since the multipliers of the hyperbolic 
fixed point of the Suris map are v and l/v. 

The series T{u) is rather intriguing. It is an analytic function of u on the 
interval {0 < u < 1}, and approaches zero rapidly as v increases. In fact, 
we will see in section 5 that r(z^) is strictly positive, but approaches zero 
exponentially fast: 



as 7/ — > 1~. This, of course, corresponds to J — > 0+. The analysis uses a 
remarkable relation between the series for F and elliptic functions, and some 
formulas found in one of the notebooks of Ramanujan (cf. [2]). 




Transport in a Perturbed Integrable Standard Map 



8 



2 Melnikov formula for twist maps 

In this section we review derivation of the Melnikov formula for twist maps 
[8]. We begin with a Lagrangian generating function S{6,6') that sat- 
isfies the twist condition (4) and has zero net flux, (5). It gives a map of 
the annulus implicitly through (3). Alternatively, the generating function 
defines the Lagrangian equations of motion through the action W defined 
on a sequence [9] = {6\ . . . , 6-^ }, by 

j 

w[e] = Y,s{0\0^^^) 

It is easy to see that an orbit of the map that begins at 6^ and ends at 0^ 
corresponds to a critical point of the action W under variation with respect 
to the interior points, and a periodic orbit is a critical point with respect 
to all points subject to the constraint that 6^ = 6^. The corresponding 
momenta are then defined through (3) as r* = —diS{6'^ ,6''^^). Thus, for 
example, a point {6a, Ta) a fixed point of the map if and only if Oa is a 
critical point of 3(6,0), and is defined through (3). 

We now describe how a modification of the Melnikov method can be ap- 
plied to predict the transversal intersection between the stable and unstable 
manifolds of two different periodic or fixed points of a twist map. This 
method is based on the variational approach of Aubry [1] and Mather [15] 
and can be applied in general to any twist map that has a saddle connection, 
in particular to any integrable twist map. We begin with a map /o gener- 
ated by 5o. Suppose that /o has two hyperbolic fixed points Za = {6a, ra) 
and = {6h,rh), and there is a saddle connection defined by the graph of 
a function x{6) on the interval U = {6a, 6t,) between these points. A diffeo- 
morphism h lU is induced by the restriction of the map to the saddle 
connection: 

M6,x{0)) = {h{6),x{h{6))) 
Let P be a function with zero net flux, (5). Then the function 

S,{6,6') = So{6,6') + eP{6,6') 

generates a twist map for small enough e. Since hyperbolic points are 
nondegenerate critical points of the action [10], the perturbed map will have 
nearby hyperbolic flxed points for small enough e. A simple case occurs 
when 6a is a critical point of P{6^6) as well as of Sq{6^6) since it is then 
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a critical point of S^{0,0) as well. Thus the fixed points of will have 
unchanged configurations, but their momenta will be modified according to 
(3). 

It is well known that there are useful relations between the action of 
orbits and areas of regions for twist maps [14]. We will use one such relation 
to obtain the Melnikov formula: a relation, between the graph x ^"^^ the 
action of orbits on the stable and unstable manifolds of a hyperbolic fixed 
point [23]. Let 6^ be a point on the unstable manifold of a fixed point Za 
that is close enough to Za so that the segment of W^{za) containing 6^ is 
given by a graph x(^)- Let < 0, be the preorbit of this point. Defining 
the backward action difference as 

t——oo 

then the unstable manifold is defined by the graph of the function 



39 



A corresponding formula for the forward action of an orbit on the stable 
manifold yields a formula for the graph of an initial segment of the stable 
manifold, x*, of Zf,: 

Aw^{0')= -Er=o[s{d\d'+')-s{e,,e,)] , 

x'ie) = ^ . (10) 

The difference between these two actions leads to the Melnikov-like for- 
mula for the transversal intersection of these manifolds. We summarize the 
results as a theorem: 

Theorem 2 Let Sq be the generating function for a twist map fo that has 
two hyperbolic fixed points Za and z^ with a saddle connection given by r = 
x{d) for 6 E U = {6a, Of,)- Let fo induce a diffeomorphism h(9) on the 
connection. Then if S(, = Sq + eP generates the twist map f^, such that the 
perturbation P has the following properties 

a) P(^„e«) = P(^6,^6) = 
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Then for e > small enough: 

i) the perturbed map has two hyperbolic fixed points continued from Za and 

zb; 

ii) the Melnikov series 

oo 

m= y: P{h\e),h'+\e)) (11) 

i=— oo 

converges absolutely and uniformly to a function on lA; and 

iii) if L has a nondegenerate critical point onU, then the unstable and stable 

manifolds of the two fixed points intersect transver sally. 
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3 Transport 

The theory of transport in two dimensional maps is based on a partition 
of phase space into regions between which transport is restricted by partial 
barriers of some sort. One of the simplest such partitions is to define a 
resonance zone associated with a saddle fixed point and its homoclinic tangle 
(or the heteroclinic tangle between two such points) [14, 20, 4]. We first 
recall a few definitions. 

An initial segment of a (un)stable manifold of a saddle z is defined as a 
segment starting at z and continuing to some endpoint, say p. A resonance 
zone for 2; is a compact region bounded by alternating initial segments of 
stable and unstable manifolds that intersect only at their endpoints. Each 
such intersection, apart from that at z defines a principal homoclinic point 
p. Even with the choice of a particular homoclinic orbit there are an infinity 
of choices for a resonance zone, since each point on the orbit of p 

is a principal homoclinic point; however, for an area preserving map these 
regions all have the same area. To be physically meaningful, the point p 
should be chosen so that the resonance zone is not too distorted; often 
symmetry can be used to give an appropriate choice. 

We will see that, for a perturbation of the Suris map, there are two 
saddles Za and Zb that have heteroclinic orbits going both from z^ Zb and 
from Zb ^ Za, see figure 2. The resonance zone is made up from four initial 
segments, first a segment of W"^{za) to p that joins a segment of W^{zb) at 
p forming the upper part of the boundary, and then a segment of W^{zb) to 
a lower homoclinic point, where it intersects a segment of W^{za) forming 
the lower piece of the boundary. 

Escape from a resonance will be slow if the unstable and stable manifolds 
are nearly coincident. Recall that the exit set is the set that leaves the reso- 
nance upon one iteration, and the incoming set is the set that enters the res- 
onance upon one iteration. These are easily obtained by taking the preimage 
of the boundary-all points on the unstable segments shrink towards their 
respective saddles, and the points on the stable segments lengthen. The 
preimage of a resonance is also a resonance, but the principal homoclinic 
points p switch to f^^{p)- For a principal point p the segment of be- 
tween p and f~^{p) together with that of connecting these points bound 
a turnstile that is the union of the exit set and the incoming set. 

For an area preserving map, the exit and incoming sets must have equal 
areas. In the simplest case there is one additional principal homoclinic point 
on the segment of the stable and unstable manifolds between f^^{p) and p, 
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Figure 2: Resonance for the perturbed Suris map with S = 0.5 and e = 0.05. 
A pair of principal homoclinic points are labeled p and q, and the exit and 
incoming sets are labeled £ and the I. 
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call it q. In this case the turnstile contains two lobes, the exit lobe and 
the incoming lobe. Specializing to the case of a twist map with generating 
function S, and denote the two principal homoclinic points by p = {rp,6p) 
and q = {vq^Oq). In this case, the areas of the exit and entrance lobes are 
determined by the difference in action between the orbits of p and of q [13] 

oo 

^= E s{e\,e\+') - sie^el^') . (12) 

i=— 00 

The area A is the signed area below the segment of between q and p 
and above that of W^. 

Now, the action is stationary on an orbit, and when there is a saddle 
connection, the actions of p and q are equal. Thus it is easy to see that the 
area of the exit set for a perturbation of a map with a saddle connection is 
determined, to lowest order, by the difference between two critical values of 
the Melnikov function L [11]. We summarize this result as a theorem. 

Theorem 3 Let Sq, P, S^, /q, /e, U, h and L as in section 2. Assume that 
Op and 6q are two points in U such that 

a) L'iOq) = L'iOp) = 

b) L"(ep) > 0, L"(eq) < 

c) L'{9) y^O for6q<6 <9p 

Then the heteroclinic points p and q continue heteroclinic points of f^, and 
the stable and unstable manifolds of f^ enclose a lobe with area 

A{e) = e{L{eq)-L{ep)) + 0{e^) 
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4 The Suris Map 

In this section we give formulas for the dynamics on the homoclinic connec- 
tion of the fixed point. Recall that the Suris map we consider is generated 
by the function (6) with the potential (7) 

Since V has minima at 6 = ^, the Suris map has hyperbohc fixed 
points at {0,r) = (—^,0). In fact, there are saddle connections between 
neighboring points and their exact formulas can be computed. Setting 
U = {—^ < < we find that the diffeomorphism on the saddle con- 
nection, hv : U ^ U m 

, ... 2 ( {u + i)tM^e) + {v-i) \ 

/ij/ ( f) = — arctan — -4— ; — , (13) 

^' TT \{v-l)i-An{^e) + {u+l) ) ^ ' 

where v = {\ - ^/5)/(l + ^/5). Then we have the following. 

Lemma 4 Let h^, be the diffeomorphism ofU given by (13). Then h^ sat- 
isfies 

a) /t^ = hyt , for all t ^Ij. 

b) K{-e) = -K-.{e). 

c) —\ is an unstable fixed point and ^ is a stable fixed point of hi,. 

d) V^{e) = h„{e) -29 + h^^{e), for eeU, for Vs given by (7). 

Proof. A direct computation proves a) and b). 

Direct computation for c) shows that the only fixed points of h^ occur 
when tan(|0) = ±1, or 9 = ±i. Furthermore, h'i,{—\) = ^ and /i'j,(^) = y. 
Since v G (0, 1) this implies that the former is unstable and the latter is 
stable. 

For d), we will first show that when 6 eU, 

. , /.X 2 / V^COSTT^ \ ^. 

0-/i40) = -arctan , (14) 

and 

g-/,;i(g) = _arctan . (15) 
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In order to prove these formulas, we notice first that ifO eU, then h{6) G U. 
Therefore -l< - h^{e) < I. On the other hand, 



This imphes (14). It is easy to see that i{ 6 & then (16) is positive and 
therefore < ^ — hi,{6) < 1. The substitution —6 ^ 9 in equation (14), 
implies that < —0 + h,^-\{6) < 1, and equation (15). We conclude that for 
all e eU, -I < 26 - h^-i{e) - K{e) < l. To finish the proof we work on 
the tangent of the second difference 



Since this last expression is — tan(|F^'(0)), this completes the proof. □ 
4.1 Saddle Connections 

With the help of lemma 4, we can give a description of the intersection of 
the saddle connection between (— 1,0) and (^,0). The saddle connection is 
shown in figure 3. 

Proposition 5 Let fg be the twist map generated by Ss- Then (— 1,0) and 
(^,0) are hyperbolic fixed points for fg and there are two saddle connections 
between them given by the graphs 




(16) 




VS cos 7r0(2\/Jsin ttO) 



1 — S sin^ it6 + S cos^ nO 
(5 siii(27r(^) 



1 + 5cos(27r6») 



and 



r = 



x-{0) 



6 



h{9). 
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r 

0.5 ■ 




1 \ 1 


\ -0.25 






0.25 >^ 














-0.5 





Figure 3: Saddle connections for the Suris map with 6 = 1/2. 
Proof. Using Lemma 4, the map can be written in the form 

fie, r) = (r + h^ie) -e + k-i ie),r + K{e) -20 + k-i (e)) 

Therefore, if r = x+(^) = ^ - K-i (0) then 

Now since 6 = — ^ is an unstable fixed point for h^,, this graph gives the right 
going saddle connection. In the same way setting r = x~(^) = ^ ~ hi,{6) 
gives 

f{e,x-{0)) = {K-i{e),K-i{0)-0) = {K-i{e),x-{K-^{e)) • 

Since the map conjugates to h^-i on this graph, this is clearly the left going 
saddle connection. □ 

4.2 Proof of the Main Theorem 



In this section we sketch the proof of the main theorem. The analysis of the 
infinite series for the Melnikov function relies on some formulae that can be 
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found in one of the notebooks of Ramanujan (cf. [2]), though we are not 
sure if he is the original author of the formulas or the first to pubhsh them. 
However, the formulas are remarkable and all that we do in this section and 
the next is a consequence of them. First we recall the theorem: 

Theorem 6 (Main Theorem) Let Vg be given by (7) and Sg be the Suris 
generating function (6). Let 

Sg^e = Sg + e cos^ ttO 

be a perturbed generating function with corresponding twist map fg^^. 

Then for < ^ < 1 there are at least two distinct transversal hetero- 
clinic orbits connecting (—1/2,0) and (1/2,0). Furthermore, the stable and 
unstable manifolds of these fixed points enclose a lobe with area A{5, e) given 
by 

where v = {1 - ^/S)/{^ + ^/S). 

Proof. Let P{6) = cos^ nO. It is clear that P satisfies the conditions of 
Theorem 2. Let 

oo 

L{e) = nh\0)) . 

t=—oo 

According to theorem 2 a sufficient condition for transversal intersection of 
the perturbed manifolds is that L has a nondegenerate critical point on the 
interval U = {—^ < 9 < The graph of L over U is shown in figure 4. 

To proceed, we define L{z) by L{6) = L(tan(|0)); therefore L{z) = 
EZ-ooMz) where 

. . 4;.^^(1 - z' f 

((l-;2)2 + I.2t(l+^)2)2 ■ 

This implies that is a non degenerate critical point of L whenever tan(|0) 
is a nondegenerate critical point of L. We are going to show that 6q = 
is a local maximum of L and 9p = ^ arctan((l — y^)/(l + is a local 
minimum. 

First, we write L as 

oo 

L{z) = aoiz)+Y,Mz) + a-t{z)} . (17) 



2.8852 



0.2 



0.4 



Figure 4: Melnikov function, L for 5 = 1/2. 

Now, since at{—z) = a^t{z), this implies that L is even, and so L'{0) = 0. 
Differentiation gives 

h (1 + ^''f 

Let 

then ^"(O) = — 8ro(^'^). Expanding the denominator of Fq and rearranging 
the sums gives 



Now let K[x) be the normalized complete elliptic integral of the first kind 
K{x) = 1 -7=^ ■ (18) 

Jo ^yl-XSm^T 

and define an increasing diffeomorphism of the interval (0, 1) 



^(-)=exp(-7r^^g^') . (19) 
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Then, a Ramanujan formula [2, III. 17, entry 14 (v)] implies that 

and therefore -^"(O) < 0, for all < < 1. This shows that 0g = is a 
nondegenerate local maximum of L, for all < < 1. 
To find the local minimum, we rewrite L as 

oo 

L{z)=Y,{at+,{z)+a^t{z)} . (20) 
Now we find that L'((l - + = and 

L" 1 - iA / 1 + = > -5^ . 

Notice that Z"((l-^)/(l + ^)) = To(i^'^)-To(i^). After the change of co- 
ordinates u = H{x), we get from formulas (III.17.14(v)) and (III.17.14(ix)) 
in [2] that 

ro{H{xf)-ro{H{x)) = K{x)Hl-x)-K{x)Hl-xf 

= K{x)^{l-x)x , 

and therefore L"((l - v^)/(l + v^)) > 0, for all < < 1. This shows 
that Op is a nondegenerate local minimum of L, for a\\ < v < 1. 
Using Theorem 3, we conclude that a lobe of area 

A{6,e) = e{L{eg)-L{ep)) + 0{e^) 

is enclosed by the stable and unstable manifolds, where 6p and 6g are given 
above. Finally, we use (17,20) to obtain 

oo 2t 



L(0,)=L(O) = l + 85:^^^^ 



and 



t=0 
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Therefore, 

oo / oo \ 

t=i \fc=i / 

oo / oo \ 

t=i \k=i / 

where r(zv) was defined in (9). This concludes the proof. □ 

The explicit formula for the first order approximation of the area is 
compared with numerical computations in Sec. 7. We conclude with a 
corollary that will be useful in the next section. 

Corollary 7 Let A{5, e) be the area of the lobe that was described before. 
Let K{x) and H{x) as in the proof of the Main Theorem. 
Let G{x) = ((1 - H{x))l{l + H{x)f. Then 

A{G{x),e) \ 
hm = KixYil - x) 

e->-0 e 

Proof, lis = G{x) then v = (1 - v^)/(l + v^) = H{x). Using [2, formula 
(III.17.14.(i))] we find that 

T{H{x)) = K{xf{l-x). (21) 

□ 
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5 Exponentially small behavior 

In this section we investigate the asymptotics of the lobe area, (8), as ^ — )■ 0. 
We will see that the area is exponentially small in S~^^'^. The analysis in 
this section is also based on the Ramanujan formalae used in the previous 
section. We summarize with a lemma. 

Lemma 8 Let be defined by (9), then 

Proof. We use (21) of corollary 7, based on rewriting F in terms of the 
elliptic integral (18) and the diffeomorphism H, (19). Our elliptic function is 
normalized so that K{Q) = 1; it has the asymptotic form [3, formula 112.01] 

1 1 

K{x) -\og{- )asi^^-l-. (23) 

TT 1 — X 



Thus H has the limits 



lim H{x) = , lim H{x) = 1 



Further, (23) gives 



H{x) ~ — as x ^- 0+ ^ H-^{^) ~ 16z^ as ^- 0+ . (24) 

The definition of H implies that logi!f(a;) logi!f(l — a;) = tt^ and therefore 
that 

and so 

1 - a; = 1 - i?"^(i^) = H-^ (exp(-7rVlog(l/i^))) . 
Combining this with (24) gives 

1 - a; ~ 16exp(-7rVlog(l/i^)) as ^- 1" . (25) 

Finally, using (25) and (23) yields 

K{x) = K{H-^{u)) ^aszv^-l". (26) 

log(^) 
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Putting equations (25) and (26) into (21) gives the promised result. □ 

The asymptotic result for ^ 1^ can be easily converted into one for 
6 ^ using 

1 1 

iog(TH ^ 

Combining this with lemma 8 gives 
Corollary 9 

^(5,6)~^exp(^^j+0(e^), (27) 

as S ^ 0. 
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6 Anti-Integrable Limit 

While it doesn't fit in with the rest of our analysis, we present in this section 
a large e expression for the lobe area. Our main purpose is to have an 
expression to compare with the numerical results in the next section. This 
expression is easy to obtain using the "anti-integrable limit" [12]. This limit 
is obtained by scaling the action by e~^, to get 

s = sie,e')/e = p{e) + e-^Ss , 

and then setting e^^ = 0. The point is that when e ^ 1 > ^ the points 
on the two heteroclinic orbits are all found in a neighborhood of the critical 
points of the potential P = cos^lnO), i.e., at 6 = m/2. In the anti-integrable 
limit, an orbit consists solely of a sequence of configuration points sitting 
at these critical points. Thus the two heteroclinic orbits are given by the 
sequences 

\"p} {■■■' ~2' ~2' 2' 2' ' ^ ii ' ~ 2 ' ~ 2 ' 2 ' 2 ' 

The action is stationary on an orbit, thus to first order in the small parameter 
e^^, the change in the orbit with e can be ignored in the action difference 
to give 

-A= Y: 0(6-2). 

t=— oo 

For the Suris map, this yields 

A = e-l-^ [dilog(l + 5)- dilog(l - 5)] + O(e-i) (28) 
We compare this result with the numerical calculations in section 7. 
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7 Numerical comparison 

In this section, we compare the theoretical results with numerical computa- 
tions of the lobe area. The task is to find the actions of the two homoclinic 
points p and q. For this task we use the symmetry of the Suris map. The 
action difference between the orbits of q and p gives the lobe area, (12). 

A map is reversible if it is conjugate to its inverse by an involution R: 
B? = I and RfR = f^^. We call R a reversor for /. Providing V is an even 
function, a reversor for the standard map is given by 

R : {6,r) ^ {-e,r +V'{e)) . (29) 

Note that if R is a reversor, then is as well. Fixed sets, Fix(i?) = {z : 
Rz = 2} of a reversor are important because 

Lemma 10 Let Za be a saddle fixed point of f . Suppose there is a point 
q G W^{za) that is fixed under R, then the orbit of q is heteroclinic from Za 
to RZa- 

Proof. Denote the orbit of g by g*, with q = q^ . By assumption lim g~* = 

>-oo 

Za- Using the reversor gives g* = = f^Rq^ = Rf^^q^ = Rq^^, thus 
lim g* = Rza- □ 

Thus to find a heteroclinic orbit from Za to Rza it is sufficient to search 
for points on the unstable manifold that intersect the fixed set of R. Such 
an orbit is called a symmetric heteroclinic orbit. Note that if g" is fixed 
under R, then g* G Fix(/^'i?) since 

/2*i?g* = f'Rf-W = f'Rq' = fV = q' ■ 

This implies that the symmetric heteroclinic orbits divide naturally into 
pairs, one has a point on Fix(i?), which we will call q and the other on 
Fix(/i2), called p. The fixed sets of the standard symmetry, (29), are 

Fix(i?) = {{e, r):e = 0}, Fix(/i?) = {{6, r) : r = 20} , 

see figure 5. To find these orbits numerically we move along W^lza) to the 
first points that intersect Fix(i?) and Fix(/ii!), respectively. For example, 
to find q G Fixi?, let = {6^,r'^) = Za + TjE'^ where E'^ is the unstable 
eigenvector, and r/ is a small parameter to be chosen below. Let + 1 be 
the first time for which the iterate of is beyond Fix(i?). Now choose a 
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Fix(R) 



Fix(fR) 




'a 



Fix(f2R) 



Figure 5: Finding some heteroclinic orbits using symmetry. Shown are the 
fixed sets Fix(/*i?), for i = 0, 1, 2 and a rough example of initial guess zq, 
leading to = 1. 

point z{6) on the line from to z^ = fz^^ parameterized by the initial angle 
6. By construction we are guaranteed that the function Z{0) = ■Kif'^''{z{6)) 
has a zero for 6 G [0°,6^]. 

We use a root finding method (Brent's method) to determine this zero 
to some precision, say p. The choice of precision influences the original 
value for rj, as well as the number of iterates until a crossing. Assuming 
W is smooth, the point z^ will be 0{rf') away from W^. After tc iterates, 
however, this error will decrease by the factor A~*'= where A is roughly the 
unstable multiplier of the fixed point. There is no sense in having this error 
smaller or larger than the precision of our root finder, so we set p ~ X^^^rf. 
On the other hand, since we start a distance -q from the fixed point, and 
wish to go a distance 0(1) to find the first crossing of the symmetry line, 
we have 77 A*"^ = 0(1). Thus, it is appropriate to set 



To find the second homoclinic point, p G Fix(/i?), we repeat the above 
analysis, using crossing of Fix(/i?) to determine tc, etc. The lobe area is 
given by the difference in action between these two orbits, from (12). 

For our computations, using IEEE double precision arithmetic, we set 




(30) 



Transpc 




Figure 6: Log-Linear plot of A/e as a function of S for various values of e 
compared with the theoretical expression (8) (thick line) and the small 6 
expression (22) (thin line) 

p = 10~^^. These computations give apparently accurate results providing 

A » 10-14. 

In figure 6 we show a comparison of the result of (8) with the numerical 
results on a log scale. The analytical result agrees well with the numerical 
results when e = 10~^. We show the same data on a linear scale in figure 7. 
Remarkably, the asymptotic formula (22) agrees with the e = 10^^ compu- 
tation within 1% up to 5 = 0.8. Note, however, that (27) provides a poor 
comparison with the numerical results since even when 6 = 0.1, log(l/z^) 
differs from 2\/S by almost 4%. In figure 7, the anti-integrable results eval- 
uated at e = 1 are also shown. We are unable to obtain numerical results 
for such a large e, as the multiplier of the fixed points is too large. 




Figure 7: Linear plot of ^/e as a function of 5 for various values of e. 
solid curve is the the anti-integrable limit (28) for e = 1. 
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8 Conclusion 

The perturbed Suris map studied here depends on two parameters, the Suris 
parameter S and the perturbation strength e. We obtained the lobe area for 
the fixed point resonance of this map for small e to 0(e) and for large 
e to 0{e^^). In the small epsilon case, we showed that the lobe area is 
exponentially small in as 5 — > 0. This result agrees remarkably well 

with the numerical calculation of lobe area even up to (5 = 0.8. 

Results showing that lobe areas are exponentially small have been ob- 
tained by now by many authors; for our purposes the most interesting results 
are those of Lazutkin for the standard map (2) [6, 5]. He shows that 

Ar-. (l + 0{lnX'^)) as k ^ . (31) 

Here InA ~ is the exponent of the hyperbolic fixed point. The constant 
ujQ ~ 1118.83 can be computed numerically to arbitrary accuracy using the 
"semi-standard" map. The exponent in the Suris map expression (22) is 
of identical form, since ln(l/zv) ~ 2\/6 is the exponent of the hyperbolic 
point in the Suris map. Furthermore, the integrable Suris map limits to the 
standard map to lowest order in S if we set k = 45, and we see that the 
exponents are the same in this limit as well. 

It is intriguing to speculate that our results could give those of Lazutkin 
in some limit, however, this is not the case, since our expression captures 
only the 0(e) term in A and neglects any terms exponentially small in e. 
If we assume both S and e are small then our map is, to lowest order, the 
standard map with parameter 

k = 4:5 + 2n'^e 

In this case the multiplicative coefficient of A would be the same in (27) as 
in (31) if we were to set 

e = !^ 0.457015 

However, the exponent in our expression is no longer the correct one, because 
we are missing exponentially small terms. So it is clear that the standard 
map is a harder problem to study than the one studied here. 

A similar Melnikov analysis is possible for other standard maps, and in 
general, for any higher dimensional twist map that has a saddle connection of 
the type described in this paper (see [7]). In [16], we proved that there exists 
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a large class of standard maps with saddle connections. In addition, other 
authors [16], [19], [22] have found examples of twist maps that are integrable. 
So, in principle, it is possible to apply our methods and formula (11) to 
these maps. 

The study of perturbations of twist maps with saddle connections in 
higher dimensions is important because it could lead to the development 
of a higher dimensional theory of transport. The question of transport in 
higher dimensions remains, and a good definition is still needed. 
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